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Abstract 

The diffusion coefficient of an inclusion in a liquid membrane is investigated by taking into 
account the interaction between membranes and bulk solvents of arbitrary thickness. As illustrative 
examples, the diffusion coefficients of two types of inclusions - a circular domain composed of fluid 
with the same viscosity as the host membrane and that of a polymer chain embedded in the 
membrane are studied. The diffusion coefficients are expressed in terms of the hydrodynamic 
screening lengths which vary according to the solvent thickness. When the membrane fluid is 
dragged by the solvent of finite thickness, via stick boundary conditions, multiple hydrodynamic 
screening lengths together with the weight factors to the diffusion coefficients are obtained from 
the characteristic equation. The condition for which the diffusion coefficients can be approximated 
by the expression including only a single hydrodynamic screening length are also shown. 
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I. INTRODUCTION 



Recent advances in experimental techniques have made the direct observation of the 
Brownian motion of /im sized objects in membranes using microscopy and imaging a routine 



process 



-7|. As a result, diffusion coefficients can be measured accurately and it is now 
possible to address the issue of the differences between Brownian motion of macromolecules 
embedded in membranes in various environments. 

Vesicles with sizes of the order of 10 /im are frequently used in experiments while the 
typical distance of a supported membrane from the substrate is of the order of 20 A 0, f|. 
Obviously, in both these general cases the coupling of membrane with its environment are 
very different and hence it will influence the Brownian motion of inclusions. In this paper, 
we investigate the influence of solvent environments on diffusion of an inclusion embedded 
in a membrane. In the biological context, there are many examples of membranes coming 
in contact with a solvent of various depth such as in tissues. 

Biological membranes can be regarded as two-dimensional (2D) viscous fluids. An impor- 



tant feature of membranes as a transport media is that they are not purely isolated [12Hl4|. 
Liquid membranes are coupled to surrounding solvents by interaction of polar head groups 
of lipid molecules with solvents; they form quasi-2D systems coupled to three-dimensional 
(3D) solvents. The coupling to the surrounding environments induces the momentum ex- 
change between the membrane and the solvents. The influence of the momentum exchange 
on the Brownian dynamics has been theoretically investigated by introducing a phenomeno- 
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hese studies have also 



22|. 



logical coupling constant or simplifying the solvents flow [8|- 
been extended to investigate the concentration fluctuations 

Despite the large number of studies, the Brownian motion of an object in liquid mem- 
branes has not yet been fully understood. In a hydrodynamic description, 2D flow in a 
bilayer membrane can be regarded as viscous and the interaction between liquid membranes 
and surrounding solvents can be taken into account by the stick boundary condition between 
them. Diffusion coefficients of macroscopic inclusions embedded in membranes were ana- 



lytically investigated for a planar membrane surrounded by solvent layers of infinite 8Mll] 



or very small thickness 12H17I]. These studies revealed that the hydrodynamic flow in a 
membrane is screened by the solvent drag force and is characterized by a hydrodynamic 
screening length. When a planar membrane is surrounded by infinite thickness of solvent, 



it is called the Saffman and Delbriick (SD) hydro dynamic screening length, and is 
given by the ratio between the 2D membrane viscosity rj and the 3D solvent viscosity rj s , 
v~ l = r]/rj s (As we shall see below, the dimension of the 2D membrane viscosity is 

that of 3D solvent viscosity times a length.) In the opposite limit of a thin solvent layer of 
the thickness h, Evans and Sackmann (ES) hydrodynamic screening length given by \fhfv is 



appropriate 



14J. In both limits, the diffusion coefficients depend logarithmically on the size 



of the inclusions as long as the size is smaller than the hydrodynamic screening length. On 
the other hand, the diffusion coefficients depend on the size of the inclusions very differently 
when the size of the inclusions exceeds the hydrodynamic screening length. These studies 
naturally lead to the interest in the hydrodynamic screening length and its influence on the 
diffusion coefficients when the solvent layer has a finite thickness. 

The solvent flow can be varied by changing the solvent thickness. The flow of solvents 
influences the membrane flow through the stick boundary condition imposed between the 
membrane and the solvents. As a result, the diffusion coefficients depend on the solvent 
thickness. The influence of the finite solvent thickness has been recently studied for dif- 



fusion of a disk 



20(, concentration 



2lN23|. correlated diffusion 



24- 



27j], and 



fuctuations 

polymer diffusion in a membrane [28[. Diffusion coefficients of other types of inclusions on 
membranes [l^, 0-32] or on Langmuir monolayers 33 1 have also been theoretically calcu- 
lated. However, the investigation on the influence of finite thickness of solvent was limited 
to numerical evaluation of the diffusion coefficients, where the dependence of the hydrody- 



namic screening length on the solvent thickness was not completely elucidated |2CH28| . In 
this paper, the relation between the diffusion coefficients and the hydrodynamic screening 
lengths are throughly investigated for an arbitrary thickness of the solvent layers on the 
basis of the analytical expression on the hydrodynamic screening lengths. 

The relation between the diffusion coefficients and the hydrodynamic screening lengths 
can be shown in a straight-forward manner for a polymer embedded in a membrane by the 



Zimm mode 
in 2D 
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where the equilibrium average of the hydrodynamic interactions is performed 



3JJ . The multiple hydrodynamic screening lengths are then found for the finite 
solvent thickness. The diffusion coefficients are expressed by the weighted sum; each term 
in the sum is a product of the weight factor and the function of the dimensionless size of the 
polymer normalized by each hydrodynamic screening length. On the basis of the analytical 
expression, the condition that the diffusion coefficient is approximately represented solely 



by the ES hydro dynamic screening length can be discussed in detail. We show that the 
diffusion coefficient cannot be approximated by the ES hydrodynamic screening length when 
both u~ l = t]/t] s and the size of the macromolecule are smaller than the solvent thickness. 

Essentially the same relation between the diffusion coefficients and the hydrodynamic 
screening lengths is obtained for diffusion of a circular liquid domain with the same viscosity 
as that of the host membrane. The diffusion coefficient of a circular liquid domain embedded 
in a membrane has been studied in relation to recently proposed raft model, where rafts are 



34-38 



It is believed 



formed by sphingomyelin and cholesterol rich liquid domains 
that rafts undergo lateral Brownian motion within a bilayer membrane and act as platforms 
for protein association and signaling 35J. Previously, the diffusion coefficient of a circular 
liquid domain of arbitrary size was derived in the limit of infinite depth of solvent layer 
or the limit of small depth of solvent layer (3, H, 39, 40]. In this paper, the results are 



generalized for the arbitrary thickness of solvent layers. The diffusion coefficient is obtained 
as a simple integral which can be expressed again as the sum of the terms given by functions 
of the same hydrodynamic screening lengths multiplied by the same weight factors as those 
for the polymer diffusion coefficients. 

In Sec. [Til the membrane hydrodynamics is reviewed. The diffusion coefficient of a poly- 
mer embedded in a membrane is obtained in Sec. IIHI In Sec. IIVI the relation between 
hydrodynamic screening length and the solvent thickness is discussed. The diffusion coef- 
ficient of a liquid domain in a membrane is obtained in Sec. |V] Finally, the last section is 
devoted to conclusions. 



II. HYDRODYNAMIC FLOW IN A MEMBRANE AND SOLVENT 

As shown in Fig. [TJ we consider the situation where the liquid membrane is supported by 
a bulk solvent on the solid substrate. The situation where the membrane is also supported by 
a solvent from above will be considered in Sec. VI. We denote the 2D flow in the membrane 
by v(r) where r = (x, y) represents a position within the plane of the membrane. The 
membrane is regarded to be incompressible, 

V-v(r) = 0. (1) 



4 




FIG. 1. Schematic picture showing a planar liquid membrane with 2D viscosity r] located at z = 0. 
It is supported on a solvent with 3D viscosity r? s . A substrate is located at z = — h bounding the 
solvent. 

Here V is a differential operator in the 2D Euclidean space. The viscous flow in the mem- 
brane can be expressed by the Stokes equation in 2D, 

7?V 2 v(r) - Vp(r) + f s (r) = 0, (2) 

where r\ is the 2D membrane viscosity, p(r) the in-plane pressure, and f s (r) the in-plane 
force exerted on the membrane from the solvent. The last quantity can be obtained when 
the solvent fluid velocities are determined. The stress tensor of the liquid membrane is given 
by 

a a p(r) = -p{v)5 a p + 7][d a vp(r) + dpv a (r)], (3) 

where 5 a p is the Kronecker delta, and a, (3 are x, y. Then Eq. (T2J) can be represented in 
terms of the stress tensor as, 

div CT + / s (r) = 0, (4) 

where (divcr)^ = ^\ da a p/dxp. 

As shown in Fig. [TJ the membrane is located in the plane at z = 0. The solvent velocities 
v( 3 )(r, z), satisfy the incompressibility condition 

V-v^(r,z) = 0, (5) 
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where V represents a differential operator in the 3D Euclidean space. We denote the 3D 
viscosity of the solvent as r) s , and the solvent flow also obeys the 3D Stokes equation, 



r) s V V 3) (r, z) - Vp (3) (r, z) = 0, (6) 

where p^(r,z) represents the pressure of the solvent. The solvent is supported on the 
substrate which is located at z — — h. The no-slip boundary condition is imposed at z = 
— h as well as between the membrane flow and the solvent flow. Through this boundary 
condition, the surrounding solvent exerts a drag force on the liquid membrane. 

The drag force in Eq. (J2J) can be expressed as f s = — (I — e z e z ) • <j^(r, 0) • e z , where e z is 
the unit vector along the z-axis. The tensorial component of I is given by Sij, and I — e z e z 
denotes the projection to the in-plane space. The stress tensor of solvent is given by 

4°(r, z) = -p^(r, z)5 ij + r/ s [«9^ {3) (r, z) + d jV f\v, z)}, (7) 

where i, j denote x, y, z. 

Using the stick boundary conditions at z = and z = — h, we solve the hydrodynamic 



equations 



be 
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rom Eq. (|5]) to Eq. ((6]) to obtain f s . In the Fourier space, f s is calculated to 

411 

f s [k] = -?7 s fccoth(A;/i)v[k], (8) 



where k = (k x , k y ) and k = |k|. The real space velocity field of the membrane flow v(r) can 
be expressed as 

/d 2 k 
— ^v[k]exp(ik-r). (9) 



The Fourier space mobility tensor G[k] associated with the velocity field is given by 



22 



33 



Ml 

k Q kp 



Ga ^ ] ~ r 1 [k 2 + vkcoth{kh)} \ 5aP ~~W)- ( 10 ) 

In order to calculate diffusion coefficients, the mobility tensor in Fourier space should be 
transformed into real space. Previously, the inverse Fourier transform of the mobility tensor 
was analytically performed only in the limits of infinite or zero thicknesses of a solvent layer. 
In the next section, the inverse Fourier transformation of the mobility tensor is analytically 
performed for an arbitrary thickness of a solvent. 
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III. DIFFUSION COEFFICIENT OF A 2-DIMENSIONAL POLYMER CHAIN 



As an illustrative example for the influence of finite thickness of solvent on the diffusion 



coefficient of a macromolecule embedded in a 2D p 
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anar membrane, we consider the diffusion 



28 



31] . Previously, the influence of the 



of a polymer chain confined in the membrane 
solvent on diffusion coefficients is analytically investigated only in the limits of very thin or 
infinite thicknesses of solvent layers. In these works, the hydrodynamic screening length is 
a key quantity in characterizing the screening of the flow of membrane by the presence of 
solvent layers. The influence of finite thickness of solvent was investigated by numerically 
evaluating the inverse Fourier transform of the mobility tensor, where the hydrodynamic 
screening length was not even defined. In this section, the hydrodynamic screening lengths 
are obtained from an analytical equation for arbitrary thickness of solvent layer. 

The conformation of a 2D polymer chain embedded in a 2D membrane is represented by 
N beads with position vectors, {R n } = (Ri, . . . , Rat), under the potential energy, 



U 



k B T 
b 2 



N 



^(Rn — Rn- 



(11) 



n=2 



where b is the Kuhn length 



421 ] . The mobility tensor associated with the beads is given by 



the inverse Fourier transform of Eq. (TTUT) as 

d 2 k 



G a /3 (Rn — R m ) 



(2vr) : 



G a/3 [k] exp [ik • (R„ - R m )] . 



(12) 



Within the pre-averaging approximation 42J, the polymer diffusion coefficient is expressed 
as 



D 



poly 



N 



dn 



N 



dm 



g(n - m), 



,o N J N 

where g(n — m) is the isotropic component of mobility tensor G a p (R 71 



Rn 



28 



(13) 



. By using 



Eq. f|T2|) . two analytical expressions for the diffusion coefficients have been derived from 



Eq. ( Tl3]) in the limits of very thin or infinite thickness of solvent layers [l7j, |28( . Here, we 
investigate the diffusion coefficient by keeping the finite depth of the solvent layer without 
taking the 



relation 



42| 



imits. By expanding l/[k + is coth(kh)} in partial fractions, we note the general 

/(*) 



dA; 



k + ucoth(kh) 



°°dk^- 



(14) 



where f(k) is an arbitrary function, kj and Cj will be later given by Eqs. (|T6|) and (|T7|) . 
respectively. By introducing Eq. (JHJ) into Eqs. f lT2|) and f fT3|) . we obtain, 



g(n — m) 



— b 2 n 2 An 
4 J 



n — m 



m 



(15) 



where Ei(— z) is the exponential integral 44]. 

In the real space, the mobility tensor is expressed in terms of an infinite number of 
characteristic lengths, kJ 1 , where Kj is determined by the following characteristic equation 



COt(Kjh) = —. 

V 



(16) 



All the roots of the equation are given by k = ±Kj with j = 1, 2, • ■ ■ . The characteristic 
lengths relative to h, l/(Kjh), depend on vh given by the viscosity ratio v = r] s /r], and 
represent the screening of hydrodynamic flow in 2D membrane due to the presence of the 
solvent. The contribution of each screening length is weighted by the factor 

2v 



3 hK 2 + hu 2 + u' 
Using Eq. (115]) . the diffusion coefficient is obtained as 

k B T 1 



(17) 



D 



poly 



4-K1] e\ 



l + ^)(21ne,+ 7 ) 
-exp( S 2 )Ei(- S 2 )" 



(18) 



where 7 = 0.5772 • • • is Euler's constant. In the above, we have defined the dimensionless 
polymer size as ej = \^NbKj/2 = R g Kj, and R g = \fNb/2 is the radius of gyration for the 
2D Gaussian polymer chain. 

The limiting expression for e\ ^ 1 is 



D 



poly 



4irr] 



lnei 



2 + ! 

2 4 



(19) 



As will be discussed in the next section, the above expression is close to the exact result 
under the additional condition of vh < 1 which is needed to replace the sum in Eq. f fT8|) 
with the term related to e\. When e\ ^> 1, Eq. (|T8|) reduces to 

k B T 1 



D 



poly 



4:711] e 2 . 



2 In e x + 7 - 1 



(20) 
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FIG. 2. (Color online) The smallest positive values for the inverse of the characteristic lengths 
against the solvent thickness. Both quantities are normalized by v = rj s /rj. The red thick line 
represents the smallest positive root of the characteristic equation, Eq. (|16p . calculated numerically. 
The long dashed line represents k/v = 1/y/uh. The short dashed line represents the result of Eq. 
(|23p . The thin solid line represents the result of Eq. (|22p . 

This expression holds regardless of the value of h as long as it is finite. The sum in Eq. f fT8|) 
can be represented by the single dominant term as long as e± ^> 1. However, the additional 
condition of vh < 1 is required when e\ -C 1, about which we shall discuss in the next 
section. 



IV. HYDRODYNAMIC SCREENING LENGTH VS. SOLVENT THICKNESS 

If the approximated diffusion coefficient obtained by taking into account only the smallest 
positive value of Kj (denoted by Hi) reproduce the exact results, then k^ 1 can be regarded 
as the effective hydrodynamic screening length. 

First, we consider the value of n\ which is the inverse of the effective hydrodynamic 
screening length as long as the higher order (j > 2) terms can be ignored. We first note the 

9 



series expansion, 

„2 



n=l 



xcota: = l + 2V - % (21) 
^ x 2 - n 2 % 2 



Since the lowest order term can be estimated as l+2x 2 / (x 2 — tt 2 ) ~ x 2 / (vh), the approximate 
expression for k\ turns out to be 

1 /2 

( 3uh + 7T 2 - V /(3z//l) 2 + 2z// i 7T 2 + 7T 4 \ 

«i « ^ — 2 J . (22) 

In the limit of uh/n 2 < 1, K\ can be further approximated as 

Ki « « (l - ^) , (23) 

where k = \Jvjh is the inverse of the ES hydrodynamic screening length defined in the 
limit of h — y 0. In Fig. HI the smallest positive values for the inverse of the characteristic 
lengths are presented against the solvent layer thickness, h. By increasing the solvent layer 
thickness h, the inverse of the hydrodynamic screening length rapidly decreases as shown in 
Fig.d 

Next we consider the condition for which the diffusion coefficient can be characterized by 
a single hydrodynamic screening length as a good approximation for the exact expression 
including multiple hydrodynamic screening lengths associated with higher order Kj. Judging 
from Eq. (|T6|) and Fig. [3] (a), Kj takes discrete values which are almost equally separated. 
When K\R g is well separated from «2-Rg and the diffusion coefficient is given by the weighted 
sum of monotonically decreasing functions of HjR g multiplied by the rapidly decreasing 
weights, the sum can be well represented by the term associated with K,iR g alone. Below, 
we show that K\R g is well separated from K2R g when R g > h and the weights rapidly decay 
when h < 1/v. 

Since we have Kj w K\ + ir(j — l)/h, the hydrodynamic screening lengths are separated 
by the factor 1/h. Hence K 2 R g is well separated from KiR g when R g /h > 1. It is convenient 
to define the cut-off size R* g = h over which the expression with K\R g could be very different 
from that with K 2 R g . 

In Fig. |3] (b), the weight factors Cj are shown against Kj/u. The weight factors Cj 
in Eq. (I17j) decrease with increasing Kj. The ratio of C2/C1 is an important factor in 
estimating whether the term related to K\ is dominant over other terms. Fig. [3] (b) shows 
that the difference between C\ and C 2 increases by decreasing the thickness of solvent layer. 

10 



(a) 




FIG. 3. (Color online) 

(a) The pictorial solution of the characteristic equation, Eq. (JT6J) ; cot(x) against x and x/{yh) 
against x for vh = 1. The cross points of lines are Xj = Kjh. The smallest positive value is x\. «i 
is obtained by K\ = x\jh (b) Cj against Kj/v. Cj represents the weight associated with each 
hydrodynamic screening length, Eq. (fT7|) . Blue squares, red circles and dots represent vh = 0.1, 

vh = 1.0 and vh = 10.0, respectively. 
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Specifically, we have Kj k\ + ir(j — l)/h and C2/C1 ~ k\Jk\ ~ vh/n 2 . Hence the diffusion 
coefficient can be approximated by the expression involving k\ alone when vh/it 2 < 1. 
When the condition vh < 1 is satisfied, C2/C1 < 1 and the diffusion coefficients can be 
approximated by those obtained by Evans and Sackmann, where C\ ~ 1/(1 + vh/2) ~ 1. It 
is then convenient to define the critical thickness of solvent h* = 1/v. If the solvent depth 
exceeds h*, the weight of C2 is not much different from that of C\. It should be noticed, 
however, that the precise estimation of the contribution from the higher modes requires the 
whole expression besides the weights. 

The expression of the diffusion coefficient depends on the kind of inclusions. As a rep- 
resentative example, we consider the diffusion coefficient of polymer to study conditions to 
use a single effective hydrodynamic screening length given by In Fig. HJ we show the 

polymer diffusion coefficients against the size of the polymer R g to study whether the poly- 
mer diffusion coefficients can be approximated by an expression without summation. When 
vh < 1, the polymer diffusion coefficients can be approximated by taking into account only 
Ki as shown in Fig. HI It is consistent with the fact that the weight C2 is smaller than C\ 
when vh < 1 since C2/C1 ~ vh/n 2 . The situation corresponds to that considered by Evans 
and Sackmann. 

When vh > 1, C± is close to C2 and the functional form of the diffusion coefficient should 
be carefully examined. When vh > 1 and R g > h holds, n±R g > 1 is satisfied. Then 
the diffusion coefficient is well approximated by Eq. ( 1201) showing 1/ (K,iR g ) 2 dependence. 
Notice that 1/ (KjR g ) 2 decays relatively fast by increasing j. When vh > 1 and the size of 
the polymer R g exceeds the solvent thickness h, the diffusion coefficient is approximated by 
the expression given in terms of n\R g alone. 

When vh > 1 and R g < h, on the other hand, a significant deviation is seen for the 
diffusion coefficients if the higher order terms are ignored, as can be seen from Fig. HI This 
deviation originates from the fact that the weak logarithmic dependence on KjR g and K2R g 
is not well separated from K\R g if R g /h < 1. Also notice that C\ is close to the other values 
of Cj if vh > 1. In such a situation, multiple hydrodynamic screening lengths should be 
taken into account. 

To summarize, we have four length scales: the critical thickness of the solvent h* = v~ l , 
the cut-off size of the polymer R* = h, the SD hydrodynamic screening length v -1 , and 
the ES hydrodynamic screening length k^ 1 (see Eq. fl22l) ). Although h* is identical to the 

12 



10 



10"' 



10" 



10 

VR 



10' 



g 



FIG. 4. (Color online)The diffusion coefficients of polymer against size for various solvent thickness. 
The solid red lines indicate the complete solution obtained from Eq. (|18p . The short dashed black 
lines indicate the approximate solution by assuming a characteristic length scale for hydrodynamic 
screening given by 1/ki. The short dashed line for z^/i = 0.1 overlaps with the exact solution. 



SD hydrodynamic screening length v , their physical meanings are different as explained 
below. 

The diffusion coefficient can be approximately expressed by using either z/" 1 or k^ 1 . For 
the thin solvent layer, h < 1/V, the weights rapidly decrease with increasing j and the 
diffusion coefficient is given by n^ 1 regardless of the macromolecule size. This is the limiting 
case considered by Evans and Sackmann. For thick solvent layers exceeding the critical 
thickness h* = 1/V, the diffusion coefficient can be approximated by the expression including 
only the single hydrodynamic screening length k± , when the size of the macromolecule is 
larger than R* = h. In this case, we have (i%2 — Ki)R g ~ R g /h > 1 and K\R g is well separated 
from K 2 -Rg- When the solvent thickness exceeds h* = 1/V and the size of macromolecules 
is smaller than R* = h, on the other hand, the diffusion coefficient shows weak logarithmic 
dependence on KjR g and multiple hydrodynamic screening lengths should be taken into 
account. The diffusion coefficient is expressed by the hydrodynamic screening length v~ l in 
the limit of h — > oo. 
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V. DIFFUSION COEFFICIENT OF A CIRCULAR DOMAIN 



In this section, we consider the diffusion coefficient of a circular liquid domain in a 
membrane (see Fig. Although the characteristic equation, Eq. f[T6"j) associated with 
the solvent and the membrane flow should not be altered by changing the inclusion from 
polymers to liquid domains, the condition that the diffusion coefficient can be approximated 
by truncating the infinite sum to a single expression depends on the size dependence of the 
diffusion coefficients. The size dependence can differ between polymers and liquid domains. 
For simplicity, we consider the case when the viscosity of the circular liquid domain is 
the same as that of the host membrane denoted by rj. Previously, simple expressions for 
the diffusion coefficients were obtained for either infinite or very thin limits of the solvent 
Here we generalize the results to arbitrary thickness of solvent layers. 



layers 
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We consider the situation for which the center of the circular object moves with the 
velocity U, and its edge is assumed to keep circular shape without any deformation. The 
velocity field inside and outside the circular domain satisfy 45] 

r,V\(v) - Vp(r) + f s (r) + F^r) ^"^ = 0, (24) 

and the incompressibility condition given by Eq. ([T]) for all r. Here f s was defined before in 
Eq. ([H]), and is the force exerted at the periphery of the circle in the direction normal 
to the circular boundary 39] . If we take the origin of the coordinates at the center of the 
circular domain and choose the x-coordinate in the direction of U, should vary according 
to the velocity U at the periphery of the circle. From the symmetry with respect to U, F^ 



can be expressed as [£ 



F w (r) =F (fo) cos#n, (25) 

where n = r/r is the outward normal unit vector at the surface of the circle of radius R, 
and 9 is the angle between U and r. 

Our task is to calculate the total force exerted on the circular domain in the steady state 



drf s + / Mtr-n, (26) 

r<R Jr=R 

where cr is the stress tensor of the liquid membrane given by Eq. ([3]) and di denotes the line 
integration. The first term represents the force exerted from the membrane flow field, and 
the second term represents the direct friction force exerted from the solvent surrounding the 
membrane to the circular domain. 
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FIG. 5. Schematic picture showing a liquid domain embedded in a planar liquid membrane located 
at z = 0. Both a liquid domain and a membrane have the same 2D viscosity rj. It is supported on 
a solvent with 3D viscosity rj s . A substrate is located at z = — h bounding the solvent in the lower 
region. 

By using Gauss's theorem, we find 

d£a-n + [ drf s 

=R Jr<R 



dr div cr + / dr f s 

r<R Jr<R 

[ drF^(r) S{r ~ R) 

Jr<R 271\R 



-F^e x r ^ cos 2 9 = -^e x , (27) 



2tt 2 

where Eq. (jl]) and Eq. are used to obtain the third equality, and e a denotes the unit 
vector along the a-direction. Equation (1271) shows that it is sufficient to calculate F^ to 
obtain the total force exerted on the circular object from the membrane flow field and the 

15 



solvent. 

The velocity field can be formally expressed as 



v a (r) = J dv'G a ,(v - rOiffV) K ^ R ' . (28) 
In real space, the mobility tensor is expressed by the Fourier transform of Eq. ( ITUl) as 

Ga0(r) = / G ^ ex P(^ k • r )- ( 29 ) 



dr' exp(zk • r)G Q/3 [k] exp(— zk • r') 
«4V>^- (30) 



Equation (1281) can be rewritten as 

d 2 k 



(2t) ! 



Let (f denote the angle between r' and U. Then we obtain 

,.8(r'-R) 



J dr' exp(-ik-r')-P^ ) (r / 



ItxR 



f 2n dtp t\ 
= / — exp [— i (k x cos (p + kySintp) R] cos tp — , (31) 

JO 27T r 

where r'/ r' is the unit orientational vector. Equation ( 130]) can be rewritten by using Eq. ( l3Tj) 
and the relation 

(dap ~ ^) r -^ = T f- k ^ (^cos^+ k y sinip) ■ (32) 

The integration with respect to ip can be performed (see Appendix for the useful relations 
to perform this integration), and the result becomes 



v r 



d k exp(zk ■ r) 



k 



(2?r) 2 r][k 2 + vk coth(fcfe)] [k 3 R 



Ji{kR)e x 

(33) 



Finally, we note k = {k cos 0, k sin 0) and r = (r cos 0, r sin #) as well as the relation 

d 2 k , r^d^ Z^dA; 

exp(zk ■ r) 



(2?r) 2 v ' 7 2tt 7 2tt 

x exp [— ikr (cos # cos + sin 9 sin 0)] . (34) 

Then the integration with respect to can be performed to obtain 

v(r)-e !°°dk MkR)Mkr) plta) m 

[) ~ X J ™2irrjkRr[k 2 + vkcoth(kh)} ' 1 ' 
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By identifying the velocity at the periphery of the domain to be U and using Eq. (j27|) . we 
obtain 

f°° Ji(kR) 2 
U = ~ F /o dk 7rrikR 2 [k 2 + vk coth(kh)]' ^ 
The friction coefficient is given by ( = — F/U. Following the Einstein relation D = k B T/(, 
we obtain the diffusion coefficient of a domain as 

JJ dom - k B i ^ ^ nvkR 2 [k 2 + ukcot h(kh)} ■ [6 ° 
This is the generalization of the result obtained by De Koker to the case of finite solvent 



depth 



39]. 



A. Limit of infinite thickness of the solvent layer 



In this limit, the diffusion coefficient of the circular object was first calculated for the 



solid circular disk by Saffman and Delbriick 



9[. For a circular liquid domain which has 



the same viscosity as the outside of the domain, the diffusion coefficient was obtained by De 
Koker 39[. The similar expression was obtained for the domain shape relaxation times 46]. 
By taking the limit of kh 3> 1 in Eq. ( fTUl) . the mobility tensor can be written as 22|, |24| 

X ( c- 



r](k 2 + uk) 



ai3 



k 2 



(38) 



If the above mobility tensor is used in Eq. (J37J), it reduces to that derived by De Koker 
In this case, the integration can be performed by using Mathematica with the use of Meijer 
G-functions [47] 



D 



knT 



dom 



1 



_r*3 2 
3/2 ^24 



[uRf 



1/2,3/2 
0,1,3/2,-1 



(39) 



This expression is useful to take the limits with respect to vR. 
In the case of vR <§C 1, the above expression reduces to 



D 



dom 



k B T 



111 \vR) 



7+ 4 



(40) 



The difference from the result by Saffman and Delbriick is the additional factor 1/4 in the 
r.h.s. of Eq. fj4TJ]) jg, 9]. This means that the diffusion coefficient of a circular domain is 
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slightly larger than that of the disk, since the flow induced inside the domain reduces the 
friction between the membrane flow and the domain periphery compared to that between the 
membrane flow and the solid edge. In the opposite limit of vR 3> 1, the diffusion coefficient 
is obtained as 



D, 



dom 



37T 2 ?7 s i?' 



(41) 



which is inversely proportional to the domain radius, R. The obtained diffusion coefficient 
is again slightly larger than that of the disk in the same limit 



10 



disk 



8 Vs R- 



(42) 



The fact that -Ddom is inversely proportional to R is consistent with the result of 2D polymer 
chain in the membrane 



B. The limit of thin solvent layer 

The diffusion in supported membranes in the vh <C 1 limit was originally considered 
by Evans and Sackmann for the solid disk immersed in the membrane 14J. The diffusion 
coefficient of a circular viscous domain embedded in the membrane was recently studied by 
us If]]. In this case, Eq. ffTO]) takes the following form 

1 



r\{k 2 



a/3 



k a 



(43) 
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where k = (v/h) 1 ^ 2 . The above mobility tensor was previously used by us 

We replace vk coth(fc/i) ~ k 2 for h — > in the integrand of Eq. (157)) . A rigorous condition 
of small h needs some care since k = [y /h) l l 2 diverges in the limit of h — > when v is finite. 
In the previous section, we have discussed the condition in detail and shown that the results 
are valid under the condition given by vh < 1. With this replacement, we obtain 

k B T 



dom 



7n](KRy 



-- h{ K R)K x {KR) 



(44) 



which coincides with our previous result 1_6| . However, it should be noted that the diffusion 
coefficient was obtained by taking into account the hydrodynamic force from the membrane 
alone in Ref. 16j. In order to compare the present result with our previous result, the direct 



friction between the solvent and the domain, nri(nR) 2 , should be added to the previous 
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result. This leads to add k B T / irr](KR) 2 to the diffusion coefficient. For comparison, we also 
write the result by Evans and Sackmann [ijj] 

-l 



D 



knT 



ES 



2 + 



AK^kR) 



ttt](kR) 2 [ kRK (kR)_ 



(45) 



where the direct friction between the solvent and the domain is added. As pointed out 



16]. This is because the fluid flow in the 



before, Eq. ( 144|) is slightly larger than Eq. (|4o| 
domain reduces the friction between the domain and the host membrane at the edge. 
In the limit of kR <C 1, the previous result is reproduced 1(| 



D 



dom 



4m] 



In 



2 

~kR 



-7 + 



In the opposite limit of kR ^> 1, the diffusion coefficient is obtained as 

k B T 



dom 



2ttt](kR) 2 ' 



(46) 



(47) 



In this limit, -Dd om decays as 1/R 2 as pointed out before [161]. 



C. Finite thickness of solvent layer 



In the case of finite h, the integration of Eq. ( 1371) can be transformed into the summations 
as employed before 



k n T 



- - Ii(KjR)Ki(KjR) 



(48) 



where Cj is the weight factor given by Eq. (I17p and Kj is determined by Eq. (I16p . 

When vh < 1, Cj decreases rapidly as j increases as already shown in the previous section. 
In this case, Eq. ( 1481) can be approximated by the lowest order expression, 



D. 



dom 



Ci 



k K T 



irrj(KiRy 



where K\ is the smallest positive value of Kj. For K\R 1, Eq. ( T49l) reduces to 



D 



dom 



c x 



k B T 
4irr] 



In 



whereas for k±R 3> 1, it reduces to 



D 



dom 



K\R 



k B T 



7+ 4 



2irr ] (K 1 R) 2 ' 



(49) 



(50) 



(51) 
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When vh > 1, Eq. (149|) approximates the exact expression, Eq. (1481) . only when the 
expression multiplied to Cj rapidly decreases with increasing j as we have already stated in 
the previous section. Since the hydrodynamic screening lengths and the weights factors are 
common, the difference between the polymer and the circular liquid domain originates from 
the non-dimensional size dependence. However, the size dependence is very similar between 
the polymer and the circular domain, i.e., weak logarithmic dependence for relatively small 
sizes and the algebraic dependence at large sizes. As a consequence, essentially the same 
results as those shown in Fig. [4] are obtained for the two cases. The condition that the 
diffusion coefficient can be approximated by the expression with a single hydrodynamic 
screening length is essentially the same for the polymer and the circular liquid domain. 
When vh > 1, Eq. f[o"Tj) is a good approximate expression of Eq. fT45|) but Eq. floUl) is not. 

Before closing the section, we compare in Fig.[6]the generalized solution of De Koker given 
by Eq. ( |37j) with the results in the two limits; the original solution of De Koker obtained in 
the limit of h — > oo and the results of Eq. (JUJ) obtained in the limit of h — >■ 0. The results 
of Eq. fl37|) with vh = 0.1 overlap with the results of Eq. f|44l) . The diffusion coefficient 
of Eq. (JUJ) is slightly larger than that of the solid disk, Eq. (j45|) . By increasing vh the 
results shift toward the original solution of De Koker. In the asymptotic limit, the diffusion 
coefficient scales with 1/R 2 in the generalized solution of De Koker, while the diffusion 
coefficient scales with 1/R in the the original solution of De Koker. In the opposite limit 
of R y , all the results of Eq. fl37j) as well as the original solution of De Koker show the 
logarithmic dependence on R as represented by Eq. (f40|) . 



VI. CONCLUSIONS 

The diffusion coefficient of an inclusion in a membrane is strongly influenced by the 
presence of solvents due to the stick boundary condition between the membrane and the 
solvent. The thickness of solvent layer is a key parameter controlling the diffusion of an 
inclusion in a membrane. In this work, the diffusion coefficient of a polymer confined in 
a membrane is obtained for arbitrary thickness of solvents. We also study the influence of 
finite thickness of solvent on the diffusion coefficient of a circular liquid domain with the 
same viscosity as that of the host membrane. Previously, the diffusion coefficient of a circular 



liquid domain was expressed by a single integral in the limit of infinite solvent thickness 
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39] 



10" 2 10" 1 10° 10 1 

VR 

FIG. 6. (Color online) The diffusion coefficients of a liquid domain against size for various solvent 
thickness. The solid red lines represent the generalized solution of De Koker given by Eq. (|37p . 
vh = 10, vh = 1, and vh = 0.1 from right to left. The result in the limit of thin solvent layer, 
Eq. (144p . overlaps with that of vh = 0.1. The red dashed-dotted line represents the asymptotic 
solution for vh = 10, Eq. f)51 1) . The red thin dashed line indicates the Evans-Sackman's expression, 
Eq. (|45p . The blue long dashed line represents the original solution of De Koker obtained by taking 
h—>oo limit in Eq. (l37|) . The blue dots denote the results of Eq. (j40j) and the blue dashed-dotted 
line indicates the asymptotic results, Eq. (jUJ). 

In this work, the integral expression is generalized to the case of finite solvent depth. The 
various analytical expressions are obtained from the integral expression. 

In general, the diffusion coefficient of inclusions decreases as the solvent thickness de- 
creases. The solvent induces drag against the membrane flow, and as a result diffusion is 
suppressed. The drag increases as the distance between the membrane and the substrate 
decreases. More precisely, the influence of solvent on the membrane flow is characterized 
by the hydrodynamic screening lengths. Multiple hydrodynamic screening lengths can be 
obtained from the characteristic equation Eq. f TIB]) . when the solvent thickness is nonzero. 
The largest hydrodynamic screening length characterizes the length scale of momentum dis- 
sipation from the membrane to the solvent. The membrane flow dissipates the momentum 
to the solvent through the stick boundary condition. For finite solvent depth, the largest hy- 
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drodynamic screening length is approximately given by ^Jhjv. The diffusion coefficients are 
expressed by the weighted sum of the functions of the non-dimensional size of the inclusion 
normalized by the hydrodynamic screening lengths. The weights are given by Eq. (fT7|) . 

By examining both the weights and the non-dimensional size dependence, we find that 
the diffusion coefficient can be approximated by the expression given by the largest hydrody- 
namic screening length alone except when both u^ 1 and the size of the inclusion are smaller 
than the solvent thickness. (The results are summarized in Table [H) When the solvent 
thickness is larger than i/" 1 , h > z/" 1 , the weights decrease slowly with increasing j. Hence 
v~ x can be viewed as the critical solvent depth below which there is a dominant mode given 
by the largest hydrodynamic screening length. If the size of the inclusion is smaller than 
the solvent thickness, K\R is close to K2R since we have K2R — K\R ~ R/h. Besides, if the 
diffusion coefficient depends weakly on KjR, the diffusion coefficient expressed by KiR alone 
cannot represent the exact summed result. In particular, the sum of logarithmic functions 
of screening lengths weighted by the similar factors cannot be represented by one of the 
terms. The situation is met when both v~ x and the size of the inclusion are smaller than 
the solvent thickness. 

The results summarized above are obtained for the finite solvent thickness. In the limit of 



h — > 00, a new length scale appears as discussed by Diamant 25[. The diffusion coefficient 
is given by the new hydrodynamic screening length v~ x . 

For relatively small inclusions, the diffusion coefficients can be approximated by the 
logarithmic function of the size normalized by ES hydrodynamic screening length \Jhjv 
when the solvent depth is small, i.e, \jv > h. In the opposite limit of h — > 00, the 
diffusion coefficient is given by the logarithmic function of the size normalized by v^ 1 . 
In the intermediate solvent depth, the diffusion coefficients are expressed by the sum of 
multiple terms and the diffusion coefficient cannot be represented by the logarithmic function 
of the largest hydrodynamic screening length. However, the diffusion coefficient may be 
approximated by a logarithmic function. This directs us to define an empirical interpolation 
of the effective hydrodynamic screening lengths by the inverse of k* — + y/vh) for any 
value of h when inclusions are small. 

The size dependence of diffusion coefficients is influenced by the solvent depth. In the case 
of a supported membrane, the typical value of h is 20 A and vh ~ 10~ 2 can be estimated by 
introducing typical values of membranes; r) s = 10~ 2 poise and r\ given by 1 poise multiplied 
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TABLE I. Summary of the size dependence of the diffusion coefficient when the solvent depth is 
finite. 



amplitude separation 


screening lengths 


solvent depth 


domain size 


size dependence 


(h < l/u) 


separation (h < R) 








complete 


a 


h < 1/ki <1/u 


R < 


ln(#sii2) 


complete 


a 


h < l/«i < 1/is 


< R 




incomplete 


incomplete 


1/v < 1/ki < h 


R<h 


b 


incomplete 


complete 


1/v < < h 


h<R 





a The summation in the expression of the diffusion coefficient can be approximated by the dominant term 

as long as amplitude separation is complete no matter about the separation of screening lengths. 
b The summation in the expression of the diffusion coefficient cannot be approximated by a single 

dominant term. Even in the regime of l/«i < R < h, 1/(kiR) 2 spatial dependence is not obtained. 

by the membrane thickness 5 x 10 -3 /im This is the case when the ES hydrodynamic 
screening length Eq. (|45j) or its modification Eq. (JHJ) are relevant. By using relatively large 
size of inclusion, R > h = 20 A , the asymptotic 1/R 2 dependence of diffusion coefficient 
can be observed. 

In the case of vesicles of 10 /im size, vh ~ 1 can be estimated by interpreting vesicle radius 
as solvent thickness 6j. Since the inclusion is normally smaller than the vesicle radius, we 
have R/h < 1. In this case, we can estimate as kR ~ uR/y/uh < \fvh ~ 1. If the vesicle 
radius is regarded as a solvent thickness, it may be difficult to observe the asymptotic 1/R 2 
dependence of the diffusion coefficient. However, this estimation is not rigorous but is done 
just for the purpose of indicating the boundary effect caused by the finite radius of a vesicle. 
The real flow inside a vesicle should be different from that in the presence of the solid 
substrate. It should be also reminded that there is an additional difficulty to differentiate 

n 

the translational diffusion of a domain from the rigid rotation of the vesicle [7[. 

For simplicity, we have considered the situation where the membrane is floated on a 
solvent layer. In general, both sides of a membrane are surrounded by solvents. We consider 
the case that the solvent layer on the membrane is covered by a substrate and is not a free 
standing film. We denote the 3D viscosity of solvent and the thickness in the upper domain 
as 77+ and h + , respectively, and those in the lower domain as r/~ and h~ , respectively. 
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Equation ({TO]) and Eq. ( )37|) are valid if we replace k + v coth(kh) with k + v + cot\i(kh + ) + 



^ coth(/c/i ), where v + = rj^/rj and v = rj s /rj 23]. Correspondingly, the characteristic 
equation becomes 

u + cot(njh + ) + v~ cot(njh~) = Kj. (52) 

The influence of solvents on both sides of the membrane can be investigated by studying 
the roots of Eq. fl52|) . By using cot a; ~ 1/x, we obtain 



when Kih + < 1 and Kxh^ < 1. 

In the simple situation where the membrane is sandwiched by the same depth of solvent 
layers h = h + = h~, both the characteristic equation, Eq. ( TT6l) and the weights given by 
Eq. (1171) still hold by using the renewed definition of v = v + + b>~ . The diffusion coefficients in 
this particular case can be obtained from those presented in this manuscript by substituting 
v = v+ + u - Q. 

If both h + and h~ are infinite, the screening length changes from k^ 1 to (z/ f + z/~)~ 1 . Our 
study on the influence of the finite solvent thickness indicates that the diffusion coefficient 
can be approximately expressed in terms of Ki except when both the SD screening length 
and the size of the inclusion are smaller than both h + and h~ . 
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APPENDIX: USEFUL RELATIONS 



We have used the relations, 



/ — exp [—i (k x cos ip + k y sin ip) a] = J (ka) (A.l) 
Jo 2tt 



24 



and 



/ — exp [—% [k x cos if + k y sin tp) a] cos 2 p 
Jo 2tt 

Ji(ka) 



(A.2) 



2tt 



2^ 



exp [— i (k x cos p + k y sin <p) a] cos sin p 



1 «9 2 



a 2 dk x dk, 



■Jo(fca) 



~^-J 2 (ka). 



(A.3) 
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